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Abstract 

The phase matrix of irregularly shaped randomly oriented large fractal 
particles is studied in the framework of the ray tracing approach for various 
refractive indices. The fractal particle model can be used for modelling local 
optical characteristics of various natural light scattering media with irregular 
particles, including oceanic water, crystalline clouds and solid aerosols. For 


nonabsorbing large particles the theory does not contain any fitting 
parameters except the refractive index of particles, which is either known 
(ice clouds) or has generally unknown distribution in a local volume of a 
scattering medium (dust aerosols and oceanic suspensions). 


1. Introduction 


Studies of photon transport in natural disperse media require 
information on their local optical characteristics (e.g. scattering 
and extinction matrices, Kokhanovsky (2001)). They can 
be easily found for the case of media with inclusions of 
a spherical shape (e.g. droplets in water clouds and liquid 
aerosols, Deirmendjian (1969)). This is not the case for most 
of the natural media, where particles of diverse (and irregular) 
shapes are often present. 

An exact solution of the Maxwell’s equations (of a Mie 
type) for an irregularly shaped particle is not possible due to 
the complexity of boundary conditions involved. The T-matrix 
approach allows for some progress (Wriedt 2002). However, 
it should be remembered that the solution of a light scattering 
problem for a single irregularly shaped particle is not very 
much helpful, if we need to find the response of an elementary 
volume of a light scattering medium with particles of irregular, 
but very diverse shapes. Clearly, the averaging procedure 
should be introduced in this case. Electromagnetic scattering 
calculations are usually very time consuming. So the averaging 
in the case of irregularly shaped particles (e.g. ice clouds 
or oceanic suspensions) cannot be performed numerically, 
considering one particle after another. Note also that the 
detailed information on the geometrical characteristics of 


irregularly shaped particles, which is needed for the averaging 
procedure, is usually not available. 

Another possibility is to use various theoretical models, 
which assume that optical properties of irregularly shaped 
particles can be modelled using combinations of particles 
of simple forms (N-particle models) or by a single particle 
of a complex shape (1-particle model). N-particle models 
are studied by Liou et al (2000) and Kahnert et al (2002) 
among others. Different 1-particle models were developed 
by Shifrin and Mikulinski (1982), Peltoniemi et al (1989), 
Macke et al (1996), and Muinonen et al (1996). For the 
sake of completeness, we note that there are also Mie-based 
models of irregularly shaped particles (Chylek et al 1976, 
Pollack and Cuzzi 1980, Drossart 1990, Grenfell and Warren 
1999). 

The accuracy of all these models cannot be checked, using 
exact theory, because such a theory is not available at the 
moment for systems of irregularly shaped particles. 

So the only possibility to find their range of validity is to 
compare them with experiments. Such a comparison is given 
here for a 1-particle model, developed by Macke et al (1996) 
for the case of randomly oriented fractal particles, which are 
much larger than the wavelength of incident light. We also 
study the dependence of the phase matrix of fractal particles 
on their refractive index. 
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2. The dependence of the phase matrix on the 
refractive index of fractal randomly oriented 
particles 


2.1. General equations 


Let us consider a 1-particle model, developed by Macke et al 
(1996). This is the so-called fractal particle model (FPM). 
In this model a system of irregularly shaped particles is 
substituted by a single fractal particle of a given complex shape. 
Then it is assumed that all essential features of light scattering 
by an ensemble of irregularly shaped particles are preserved 
in light scattering characteristics of this single particle at a 
random orientation. 

We are interested in study of the angular distribution of 
scattered light intensity and polarization characteristics for 
ensembles of highly complex particles. By definition, these 
distributions should be almost featureless due to a chaotic 
nature of scattering in this case. These smooth distributions 
also arise for a single particle, but of a complex shape. 
However, there is no a priori assurance, that the model 
gives accurate quantitative results for a given complex light 
scattering medium. Therefore, the correctness of a given 
1-particle model can be checked only by comparison with 
experiments. We make such a comparison in the next section. 

Here we present results of numerical calculations of the 
phase matrix in the framework of the FPM. The exact structure 
of the irregularly shaped fractal particle (and the code used for 
calculation) is described elsewhere (Macke et al 1996). 

The particle is constructed in the following way. The 
initial shape is a tetrahedron. A part of its triangular surface 
is replaced by a reduced version of the same tetrahedron. 
The resulting body is called the first generation of the triadic 
Koch fractal. Repetition of this procedure at the smaller 
triangles leads to higher generations. We will consider here 
the Koch fractals of the second generation. Optical properties 
of Koch fractals of the second and higher generations differ 
insignificantly (Macke et al 1996). 

The interaction of light with a large fractal particle is 
studied in the framework of the geometrical optics, using 
the Monte-Carlo ray-tracing technique. Results are averaged 
with respect to the particle random orientation. The peak in 
the forward direction is described in the framework of the 
Fraunhofer diffraction theory (Macke et al 1996). 

We also assume that all planes of a fractal particle have a 
roughness. The roughness is modelled in the following way. 
For each reflection—-refraction event the normal of the crystal 
surface is tilted randomly around its original direction. The 
azimuth tilt angle is chosen randomly with equal distribution 
from the interval [0, 27r]. The zenith tilt angle is chosen from 
the range [0, 7&/2], where & is the roughness parameter. We 
assume that € = 0.3 in this study. In principle, € € [0, 1]. 
It was found by Macke et al (1996) (see their figure 12) that 
larger values of é influence the result insignificantly. 

We assume that the wavelength is 0.54m and the 
refractive index n varies in the range 1.1-1.5, which 
corresponds to most of the natural dispersive media in the 
visible region. The length D of an initial tetrahedron’s side is 
fixed in this study and is equal to 100 um. We also introduce 
the size parameter x = kV/S ~ kD > 1, where k = 27/å, 
V is the volume of a fractal particle and S is its surface area. 
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The particle is assumed to be nonabsorbing. It means 
that its normalized phase matrix does not depend on the size 
in the geometrical optics approximation (except at the small- 
scattering region, where Fraunhofer diffraction takes place). 
This makes model not very flexible as far as real experimental 
data (which can vary for different media) are of concern. On 
the other hand, there is an advantage in this. Namely, we have 
a fixed phase matrix for a chaotic scattering and measured 
matrices for various media can be checked against this simple 
case of a single fractal particle without introducing any fitting 
parameters (if the refractive index is known and there is no 
light absorption). 

Let us describe the results obtained. For this we 
need to introduce the phase matrix P. This matrix relates 
Stokes vectors of incident Sa Gai So2, S03, So4) and scattered 
S (S1, S2, $3, S4) light with the following matrix equation: 


S=aPSp, (1) 


where @ = Osca J4nr?, Osca 18 the light scattering coefficient 
and r is the distance from a scattering point to the observation 
point. We will assume that we have an ensemble of identical, 
but randomly oriented fractal particles. Then the matrix Ê 
takes the following general form (Kokhanovsky 2001): 


(2) 


Note that in contrast to the case of spherical particles P22 # Pi, 
and P33 ~ P44. Also the parameter 


q = y Pint P33 + P34 (3) 


is not equal to | as it is the case for identical spheres. Here we 
introduced normalized matrix elements: pj; = P;j/ P11. 

Equality g = 1 means that polarized light scattering 
cannot produce unpolarized light. To illustrate the physical 
meaning of the parameter q, we consider the illumination of 
a scattering medium with the phase matrix (2) by a linearly 
polarized light with the Stokes vector 


(4) 
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which corresponds to the linear polarization at the angle 45° 
to a scattering plane. Then we have from equations (1), (2) 
and (4): 


S=a : (5) 


It holds for a completely polarized scattered light (Rozenberg 
1955): 
S = 83+ S$ +S3 (6) 


or (see equation (5)) 


Phy = Ph + Ph + Pay. (7) 
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Therefore, if equation (7), which is equivalent to g = 1, holds, 
the unpolarized light is not produced in a light scattering event. 
This is not the case for fractal particles. 

Also for spheres (both polydispersed and monodispersed 
ones) so-called cross-polarization ratio (Kokhanovsky and 
Jones 2002) 

_ l= Pp» 

1—2pi2 + pr 
is equal to zero. Note that r is equal to ratio of intensities 
registered for two experimental situations, namely for the case 
of a disperse medium placed between crossed vertical and 
horizontal polarizers, Jyy, to the case of the same medium 
placed between horizontal polarizers, [yy (Kokhanovsky and 
Jones 2002). Clearly, for a spherical particle incident vertically 
or horizontally polarized light does not change the state of 
polarization due to the scattering event. It remains either 
vertically or horizontally polarized. So there will be no light 
coming from the system of spherical particles, placed between 
crossed polarizers. This is not the case for irregularly shaped 
particles (Card and Jones 1999). 

Note, that both g and r do not depend on the size 
of particles in the framework of the geometrical optics 
approximation (Kokhanovsky and Jones 2002). However, they 
do depend on their refractive index and shape. 


r 


(8) 


2.2. The phase function 


Let us discuss results of numerical calculations now. The 
first element of the phase matrix Pı; is given in figure 1 
at various values of the refractive index n for the case of a 
fractal particle, described above. The dependence P}; on the 
scattering angle 6 is called the phase function. It describes the 
conditional probability of light scattering in a fixed direction 
and normalized as follows: 


1 Pid 
>| P11(0) sino dé = 1. 


Ih (9) 


It follows from equation (9) that for isotropic scattering, 
when scattering does not depend on a scattering direction, 
we get: Piy = 1. Fractal particles with refractive index 
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Figure 1. The phase function of fractals for various values of real 
part of the refractive index. 


n = 1.5 have the phase function, which is almost isotropic 
in the backward hemisphere. Such phase functions were 
experimentally measured for dust aerosols (von Hoyningen- 
Huene and Posse 1997, Volten 2001, Volten et al 2001). 

Phase functions in figure 1 are almost featureless and have 
smaller values for smaller n at scattering angles larger than 
x /4. For smaller angles (0 < 7/4), the opposite is true. Note 
that oscillations in figure 1 (and in figures, which follow) are 
due to a statistical noise of a Monte-Carlo code used. 

It should be emphasized that P;;(@) for systems of fractal 
particles considered do not depend on the refractive index 
of particles around 6 œ~ 45° (see figure 1). Note that it 
may explain the existence of a cross-point at this angle for 
experimentally measured oceanic phase functions (Tyler 1961, 
Zaneveld and Pack 1973). It is known that oceanic water 
contains a great portion of irregularly shaped particles, having 
various refractive indices (Shifrin 1988, Kokhanovsky 2001). 
Also oceanic phase functions are featureless as those given in 
figure | (see, e.g. the case n = 1.1). 

The case n = 1.3 corresponds to ice fractal particles in 
visible and well reproduce crystalline clouds phase functions 
(see, e.g. Macke et al (1996)). The case n = 1.5 roughly 
corresponds to dust particles (Volten 2001, Volten et al 2001). 

The asymmetry parameter, which is defined as 


1 T 
g= J) Pi; (8) cos 0 sin 8 dé, (10) 
0 


was also found, using data given in figure | and calculations 
at other values of n (see table 1). 

Note that this parameter can be presented in a following 
general form for large nonabsorbing particles (Kokhanovsky 
2001): 

1+ 80 

== 
where go is the asymmetry parameter related to the geometrical 
optics scattering process. This process leads to almost full 
isotropization of photon scattering directions at large n (see 
figure 1). Thus, we get go —> Oasn — oo. This is confirmed 
by data given in table 1. Note that for ice crystals in the visible 
region (n = 1.31) we get g = 0.74, which coincides with the 
experimental result, obtained from measurements, performed 
by Garett et al (2001) in natural crystalline clouds (see the 
average asymmetry parameter in table 4 of paper by Garett 
et al (2001)). Is it only by chance or due to the possibility of 


(1) 


Table 1. The asymmetry parameter. 


n Eo &§ 
1.0 1.0 1.0 
1.1 0.86 0.93 
1.2 0.68 0.84 
1.31 0.48 0.74 
1.4 0.39 0.69 
1.5 0.30 0.65 
1.6 0.23 0.61 
1.7 0.18 0.59 
1.8 0.15 0.57 
1.9 0.12 0.56 
2.0 0.10 0.55 
oo 0.0 0.50 
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the FPM to capture in a correct way light scattering by complex 
crystalline media, remains to be seen. 

The asymmetry parameter is around 0.93 atn = 1.1, 
which is typical for oceanic hydrosols with particles having 
almost the same refractive index close to that of water (Shifrin 
1988). Clearly, we get go > 1 as n — 1 (see table 1). 

The following parametrization holds with error less than 
5% atn < 1.9: 

go = 0.1 + 1.2exp(—Y(n — 1)”), (12) 
where Y = 2.65, no = 0.68. Equations (11) and (12) can be 
used for the modelling of the radiative transfer in ice clouds, 
using selected approximations (Kokhanovsky 2001), which do 
not require information on the phase function. 


2.3. The normalized phase matrix 


Let us consider now the normalized phase matrix elements 
Pij = Pij/Pi1. We start from the element pj2. The degree 
of polarization of initially unpolarized light after scattering 
of light by a fractal particle is equal to — pı2. The element 
Pi2 is given in figure 2(a) for various n. It is negative 
for most of the angles. This means that oscillations of an 
electric vector are predominantly in the plane perpendicular 
to the scattering plane. Generally, curves in figure 2(a) are 
similar to the case of Rayleigh scattering but with strongly 
reduced values (e.g. it is only 5% polarization at maximum 
atn = 1.5 as opposite to the case of Rayleigh scattering, 
where the maximal polarization is equal to 100% at 6 = 90°). 
Note also the shift of the maximum polarization angle Omax to 
larger angles as compared to Rayleigh scattering. Very small 
values of the degree of polarization and curves similar to that 
in figure 2(a) were found experimentally for oceanic water 
(Beardsley 1968), ice clouds (Dugin and Mirumyants 1976) 
and Sahara sand particles (Volten et al 2001). Generally, the 
degree of polarization decreases with n. Larger values of n 
lead to larger degree of randomization of photon polarization 
states due to a scattering event. 

We see that fractal particles do not produce a strong 
polarization of incident unpolarized light, which is a general 
feature of disperse media with irregularly shaped particles. 

The element pz is given in figure 2(b). It is equal to | for 
spherical particles at any 6. So the difference A = 1 — py 
can be used as a measure of the particle’s nonsphericity. We 
see that deviations from the spherical particle case (p22 = 1) 
grow with n. This result is quite understood. Namely, for 
optically soft particles there are almost no shape dependences 
of light scattering characteristics (Sid’ko et al 1990). Note also 
minima in figure 2(b) at n = 1.1 and 1.5. Such minima were 
observed by Voss and Fry (1984) and Volten (2001) for oceanic 
water and quartz aerosol particles, respectively. Generally, A 
is larger at the backward hemisphere. 

The element p33 is given in figure 2(c). We see that curves 
p33(0) for all n (except n = 1.5) are close to each other around 
scattering angle 120°, where the value of p33 is small. The 
absolute values of p33 decrease with the refractive index in the 
forward hemisphere. 

The element p34 is given in figure 2(d). It takes small 
negative values. This element describes the effectivity of the 
transformation of linear polarized light to circular polarized 
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light for a given scattering process. Such a transformation for 
fractal particles has a low probability as seen from figure 2(d). 
Generally, p34 decreases with n. 

The element p44 is given in figure 2(e). Note that the 
case n = 1.1 differs considerably from the case n > 1.2. 
The value of p44 describes the reduction of the degree 
of circular polarization for initially completely circularly 
polarized incident light beam. 

Angular dependences q(0),r(0) (see equations (3) 
and (8)) are given in figure 3 atn = 1.31 andn = 1.5—10.005, 
which roughly correspond to ice clouds and dust aerosols, 
respectively. We see that the cross-polarization r by ice 
crystals is larger than that for dust particles in the backward 
hemisphere. This is mostly due to the effect of light absorption 
by aerosols, which diminishes the influence of the particle 
shape on backscattering. 

The value of q is generally lower for ice crystals in the 
backward hemisphere. It remains smaller than 0.3 in this 
angular range. Remind that g = 1 for identical spheres. 
We see, therefore, that backscattering region is potentially 
informative for the different particulate media discrimination 
(Card and Jones 1999). 


3. The comparison with experiments 


3.1. Ice clouds 


Let us compare results obtained from the FPM with 
experimental data, starting from the case of ice clouds. Ice 
clouds are composed of particles of various shapes. Their 
phase matrices were measured, e.g. by Dugin and Mirumyants 
(1976). Results of comparison of their experiments with phase 
matrix elements at n = 1.3, which is close to the refractive 
index of ice in the visible region, where measurements 
were performed, are given in figure 4. Let us analyse them. 

We see that the measured and calculated functions pı2(0) 
have similar shapes, producing a broad negative maximum 
at side scattering angles. The absolute values, given by the 
experiment and the theory, differ, however. This comes at no 
surprise. Indeed, the degree of polarization in measurements 
is obtained as a ratio (ij — i2)/(i; + i2), where i; and iz are 
scattered intensities for two orthogonal polarizations. They 
are small and measured with estimated error up to 20%. This 
means that their differences, which are small numbers, can 
have errors exceeding 40% and more, which are not acceptable 
for any comparisons with the theory. Also crosses represent 
average over five experiments. Actual experimental data are 
scattered almost uniformly in the range pı2 € [—0.2, 0.0]. 

The accuracy of the FPM for elements p33 and p34 is quite 
high, especially taking into account that no fitting parameters 
have been used. The sign of p34 was different in the experiment 
due to a different definition of the Stokes vector p34 by Dugin 
and Mirumyants (1976). So we multiplied experimental data 
for p34 by j = —1 and presented these modified data in 
figure 4. 

Differences between calculated and measured values of 
p44 and pn are quite high. However, we note again that crosses 
(for p44) and points (for p22) present only averaged data. Real 
experimental data are scattered around these averaged curves. 
The theoretical curves both for p44 and pz are practically 
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Figure 2. The same as in figure 1 but for the element p12 (a), P22 (b), p33 (c), p34 (d), pa, (e). 


inside the experimental data scattering area, given by Dugin 
and Mirumyants (1976). Interestingly, measured values of p22 
are closer to calculated values of p44 and vice versa. We do 
not have an explanation for such a peculiarity. 

Overall, the FPM describes (at least qualitatively) main 
features of polarized light scattering by ice clouds. It can 
be used as input to the vector radiative transfer equation for 


studies of polarized radiative transfer in realistic crystalline 
cloudy media. Spherical particle models cannot be used in 
this case. Interestingly, the FPM has no fitting parameters for 
crystalline clouds in the visible region. 

Takano and Liou (1995) compared experimental results 
of Dugin and Mirumyants (1976) with theoretical ray-tracing 
calculations, assuming the following forms of particles: solid 
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Figure 3. The dependence of parameters q (a) and r (b) on the 
scattering angle for ice clouds (n = 1.3) and dust aerosols 

(n = 1.5 — 10.005), calculated in the framework of the FPM at 
D = 100 um, à = 0.5 um. 


plates and dendrites. They stated that the case of dendrites is 
closer to experimental data. Although, the theoretical data, 
given by Takano and Liou (1995) are generally closer to 
the experimental results than the FPM calculations presented 
above, their theoretical curves have special features, which are 
absent both in the experiment and in the results, obtained with 
the FPM. Thus, the FPM captures main effects of the shape 
irregularity on cloud polarization characteristics in a correct 
way as compared to modelling ice cloud particles by plates or 
dendrites. N-particle models (Liou et al 2002) can provide 
a better overall accuracy just because they have more fitting 
parameters involved. Fitting parameters are virtually absent in 
the FPM, however. This is the main advantage of the model, 
considered here. 

In conclusion, we present the comparison of the 
phase function derived in the framework of the FPM with 
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Figure 4. The comparison of calculations (lines) and measurements 
(@,+, A, E, x: see paper by Dugin and Mirumyants (1976)) of the 
normalized matrix elements. Points and crosses give measured 
values Pom and p44m, respectively. Squares and triangles 
correspond to measured values of p34 and p33, respectively. 


—— theory 
@® measurements 
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Figure 5. The phase function of an artificial crystalline cloud, 
obtained from measurements (Barkey and Liou 2001) and the FPM. 


experimental measurements, performed by Barkey and Liou 
(2001) (see figure 5). The theoretical results and experimental 
data correspond to each other quite well. The difference can 
be attributed to the fact that a lot of small crystalline particles 
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(and, possibly, spherical droplets) were present in the chamber 
during the experiment. Measurements for natural clouds give 
featureless phase function in the backward scattering region 
(Baran et al 2001), which well corresponds to results obtained 
from the FPM. 

We underline that due to the absence of light absorption the 
scattering characteristics in the visible region do not depend on 
the size of crystals in the framework of the geometrical optics 
approximation (A/D — 0). This explains why our model 
of a single particle with D = 100 um works quite well in 
the description of light scattering characteristics of crystalline 
clouds, characterized in reality by particle size (and shape) 
distributions. 


3.2. Oceanic water 


Oceanic water is another example of media with irregularly 
shaped particles (Shifrin 1988, Kokhanovsky 2001). The 
modelling is much more complex in this case. First of all, 
oceanic particles have different refractive indices (usually in 
the range 1.05—1.2), which complicates the analysis as opposite 
to the case of a well-known refractive index for ice. Also 
both large and small size particles play a role. Moreover, 
the Rayleigh scattering and a priori unknown light absorption 
by particles and water itself should be accounted for in the 
calculation of the phase matrix. For the sake of simplicity, 
we, however, neglect all these features and present results of 
comparison of averaged phase matrices of oceanic water (Voss 
and Fry 1984) with results of calculations at n = 1.1 in the 
framework of the FPM in figure 6. 

We have better overall agreement of the experiment and 
the theory at scattering angles larger than 120°. For smaller 
angles, however, the difference is considerable (especially 
for elements pı2 and p33). Voss and Fry (1984) found that 
P33 © p44 and p34 © O for oceanic water. This is consistent 
with a FPM as well (see figure 4). 

Note that small values of pı2 (similar to those, which 
follow from the FPM) are also reported in literature (Beardsley 
1968). General dependences p12(0) and p33(@), obtained from 
the fractal model at n = 1.1 and the experiment, are similar, 


== 
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Figure 6. Normalized matrix elements, obtained from 
measurements (@, A, *: Voss and Fry (1984)) and calculations at 


n = 1.1 (@—px, A—p33, X—P12). 


but actual values differ, which could be due to the influence 
of Rayleigh scattering by water and other factors, that we 
have mentioned in the beginning of this section. The FPM 
data presented in figure 6 can be used as input parameters 
for the theoretical vector radiative transfer studies in ocean. 
The advantage of this model is that it is defined by just one 
number (n). The correspondence of the model to experimental 
data can be improved, taking into account other parameters, 
which influence light scattering in ocean. This, however, is 
well outside the scope of this paper. 


3.3. Quartz particles 


We conclude the investigation of the validity of the FPM as 
compared to experiments, studying the case of yet another 
medium with highly irregular particles—dust aerosols. The 
main component of the dust aerosol over deserts is well 
represented by quartz aerosol particles. Measurements of 
phase matrices of quartz particles in the visible region were 
reported by Volten (2001). The comparison is given in figure 7, 
where we assumed that the imaginary part x of aerosol quartz 
particles is equal to 0.005 and n = 1.4. The correspondence 
of the experiment and the theory is quite good for all matrix 
elements (except p44 in the forward hemisphere). Note that the 
value of p34 is positive as compared to the case of nonabsorbing 
aerosols (compare figures 4 and 7). So the sign of the element 
p34 (if properly defined) can serve as an indication of light 
absorption inside scattering particles. 

The poor correspondence of the experiment and the theory 
for the element p44 can be explained by the difference of actual 
quartz refraction and absorption coefficients from assumed 
values. 

Phase functions obtained from the theory and experiment 
(Volten et al 2001) for quartz particles are compared in figure 8. 
Again we see an excellent agreement of the FPM to the 
experimental data. 

Volten (2001) and Volten et al (2001) gave also the average 
aerosol phase matrix, taking into account measurements 
of phase matrices by various samples, including Feldspar, 


matrix elements 
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scattering angle, degrees 


Figure 7. Normalized matrix elements, obtained from 
measurements (M, x, x, O, ®©: Volten et al (2001)) and calculated 
ones at the complex refractive index 1.4 — 10.005 and D = 100 um, 
à = 0.5 um. 


A A Kokhanovsky 


1000 


100 


10 


phase function 


0.1 


o 20 40 60 80 100 120 140 160 180 
scattering angle, degrees 


Figure 8. The phase function, normalized at the scattering angle 
30°, obtained from measurements (@: Volten et al (2001)) and 
calculated at the complex refractive index 1.4 — i0.005 and 

D = 100 um, A = 0.5 um. 


matrix elements 
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Figure 9. The same as in figure 7, but for the averaged normalized 
phase matrix of aerosol particles (Volten et al (2001)). Theoretical 
curves differ from those in figure 7 only due to the different 
statistics. We used 1000 orientations of the fractal particle for this 
figure as opposite to 10000 orientations of the particle in figure 7. 
The number of rays used was 1000 for both cases. 


Red Clay, Lokon Ash, Loess, Pinatubo Ash, Quartz and 
Sahara dust at wavelengths 441.6 and 632.8 nm, which totals 
to 14 measurements. This matrix is given in figure 9 together 
with results obtained from the FPM, which are identical to that 
given in figure 7. The agreement of this average matrix model, 
advised for use in aerosol remote sensing problems (Volten et al 
2001), with the FPM is rather good (again with the exception 
of the element p44). 

Phase functions obtained from the theory and averaged 
experimental data (Volten et al 2001) for the aerosol are 
compared in figure 10. Also we see a good agreement of the 
FPM and the experimental data. 

This suggests that the FPM may appear to be useful in the 
aerosol optics studies as well. Note that Volten et al (2001) 
compared their measurements with randomized Gaussian 
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Figure 10. The same as in figure 8, but for the averaged phase 
function of aerosol particles (Volten et al 2001). The theoretical 
curve does not differ from that in figure 8. 


particles’ model (GPM). They found that the GPM describes 
quite well their data (although theoretical calculations give 
higher values of element p44 as compared to the measurements 
(see figure 13 in paper by Volten et al (2001)). This means that 
both FPM and GPM (with a special choice of parameters) give 
similar results. Taking into account their inherent differences, 
we summarize that what matters is not a particular 1-particle 
model, but a degree of irregularity in 1-particle model chosen. 

One point should be mentioned. Aerosol particles (unlike 
ice crystals) can absorb light in the visible region. It means that 
light scattering characteristics depend on the size of particles 
even in the geometrical optics approximation (and generally 
on the particle size distribution). This may appear as a severe 
constraint in the value of the present analysis, which rely on 
calculations for a single particle with D = 100 um. However, 
in reality, what plays a major role is not a specific particle 
size distribution but rather a product 3 = x (x), where (x) 
is the average size parameter in the scattering ensemble of 
particles. This can explain a fairly good agreement we have 
obtained, comparing our model of a single particle (having 
the size parameter x = (x)) with aerosol experiments. Data 
for the value of x for the particles in the experiment are not 
available. However, we think that the value of 3 used in our 
calculations was comparable to that in the experiment, which 
explains a fairly good correspondence of the experiment and 
numerical calculations. 


4. Conclusions 


In conclusion, we have studied here the dependence of the 
phase matrix of fractal particles on their refractive index. It was 
found that this model describes fairly well the light scattering 
by systems of large irregularly shaped particles. The theory 
has no fitting parameters for nonabsorbing particles, providing 
that the refractive index of particles is known. It means that 
systems of large irregularly shaped particles (as far as their 
optical properties are concerned) differ from one another only 
by the value of the refractive index in the framework of the 
FPM. This statement is only approximately correct, because 


Optical properties of irregularly shaped particles 


different particles may have different morphology, which is 
clearly the case for dust and ice clouds. 

The model studied provides useful data for theoretical 
investigations of polarized light propagation in natural media, 
composed of irregularly shaped particles. For this one should 
average unphysical oscillations in curves presented above, 
which are due to a statistical noise of a Monte Carlo method. It 
can be done either increasing the number of rays/orientations 
or using various smoothing procedures. 

The discrepancy of the theory and experiment for the 
element p44 of quartz particles should be clarified in a future 
research. 

We hope that this paper will attract more attention to the 
FPM, especially with respect to atmospheric scattering studies. 
It is known that the Mie theory gives zero values for differences 
Pit — P22 and p44 — p33, which is not the case for irregularly 
shaped particles (e.g. for solid aerosols and ice crystals in the 
terrestrial atmosphere). Therefore, the spherical particle model 
cannot be applied to describe polarization characteristics of 
light interacting with disperse media, having a large portion of 
nonspherical particles (even for rapid estimates and qualitative 
studies). More close to reality models should be used in 
this case. One of such models is studied here, which is 
confirmed by its fairly good performance as compared to 
experimental data. 
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